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Abstract 

We present in this work a generalization of the solution of Gorenstein and Yang to 
the inconsistency problem of thermodynamics for systems of quasi-particles whose masses 
depend on both the temperature and the chemical potential. We work out several solutions 
for an interacting system of quarks and gluons and show that there is only one type of 
solution that reproduce both perturbative and lattice QCD. 
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1 Introduction 



In the past two decades, strongly coupled matter under extreme conditions has been object of 
intensive study, as it must describe from the core of some neutron stars to important features of 
the early universe and heavy ion collisions experiments (see e.g. PQE])- It is the understanding 
that hadrons are composed by asymptotically free particles that led to the idea that at suffi- 
ciently high temperature T, and/or quark chemical potential /i, the hadronic matter should be 
described by quarks and gluons degrees of freedom [3], the quark- gluon plasma (QGP). Current 
lattice calculations strongly suggests that this phase transition actually happens, in the case 
of a vanishing chemical potential, at a temperature around T c ~ 190MeV [I]. The search for a 
phase transition in the case of a non- vanishing quark chemical potential has recently advanced, 
as a result of new techniques for lattice calculations [51 El E] ■ 

Results from ultra- relativistic heavy ion collisions at RHIC [HI El [101 H] , on the other hand, 
indicate that the QGP has already been produced, and that hydrodynamics is capable to de- 
scribe the experimental data [12]. The basic input for hydrodynamics is the equation of state 
EoS, which must cover both the hadronic and the QGP sectors. The matter created in the 
heavy ion collisions has high T and small /i, therefore the study of the EoS at these conditions 
is of extreme importance for the hydrodynamics approach. The EoS of the hadronic phase is 
currently described by the hadronic resonance gas model [13J, with a good agreement to lattice 
data [SHU] - For the deconfined phase, one needs an EoS in terms of quark-gluon degrees of free- 
dom. QCD at finite T and fi is, in principle, the theoretical tool to compute the thermodynamics 
functions in this new phase. However, strict perturbation theory, which have been pushed to 
g 6 ln(l/g) [15] . is reasonable only for extremely high temperatures; at temperatures near T c it 
seems not to be applicable, and alternative treatments appear to be necessary. Moreover, the 
perturbative series seems to be weakly convergent [T5], [T6], [TT| . Specifically, it is expected that 
when T ^> T c , the plasma behaves like an ideal gas of quarks and gluons, but the perturbative 
series converges at a slow pace towards this expectation. A way to circumvent this problem 
is through the reorganization of the perturbative series. For instance, there are attempts to 
reorganize the series using the Hard-Thermal-Loop (HTL) effective action [18], in the form of 
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so-called HTL perturbation theory [121 ED], and also attempts based on the 2-loop ^-derivable 
approximation [211 [22]. The latter approach assumes a massive quasi-particle formalism [23, 24]. 

Concomitantly, phenomenological models based on quasi-particles have had success in de- 
scribing the lattice data [231 1211 [251 I2S1 [271 1211 [291 [30], over a wide range of temperatures (and 
at n — 0), in addition to being compatible with the fluid-like behavior found at RHIC [31J. 
More recently, a quasi-particle model for the quark-gluon plasma (qQGP) for a non-vanishing 
/x has also been worked out [321 [331 El], with a satisfactory description of the lattice data. 
Besides the quasi-particle approach, there are other treatments to describe lattice data near 
the transition region as, for instance, the Polyakov loop enhanced Nambu-Jona-Lasinio (PNJL) 
model [351 EH], which is able to fit lattice data quite well, even when /i ^ 0. 

The quasi-particle model of the quark-gluon plasma is a phenomenological model that as- 
sumes non-interacting massive quasi-particles which, with the aid of few parameters, like the 
thermal masses, is able to fit lattice QCD data over a wide range of temperatures [2^1 ESI I2"7| : 
Not only at extremely high temperatures, as in strict perturbation theory, or when T > 20T C , 
as in the HTL perturbation theory, but also near T c . Indeed, quasi-particle models intend to 
describe deconfined matter from T c up to T — > oo, although from the point of view of linking the 
QGP with hydrodynamics calculations, the reproduction of lattice data near Tc is somewhat 
more important. 

The problem with the quasi-particle models is that the thermodynamics relations calculated 
from them are not consistent [37J- Its origin is in the use of Standard Statistical Mechanics (SSM) 
[38J with dispersion relations for quasi-particles that resembles the dispersion relations for free 
particles. Contrary to the case of particles, quasi-particles have masses that are dependent 
on T and/or fi. The first work to systematically study the consistency problem, and solve it 
for finite T but zero chemical potential, was the one by Gorenstein and Yang [37]. In their 
approach, it is built an effective Hamiltonian which is composed of a part representing an ideal 
gas and a part representing contributions from the vacuum. This extra vacuum term leads to a 
modification in the pressure and in the internal energy of the system, but keeps the entropy and 
the number density untouched. In this way, a consistent thermodynamics for quasi-particles is 
obtained. Most of the works in the literature are based on this approach [39J. More recently, 
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several authors have also turned their attention to the thermodynamics consistency of quasi- 
particle systems. In [ID], the authors build an effective Hamiltonian, H e ff, as before but require 
that the average of dH e ff/dfi and dH e jf/dT are zero. As a result, they get expressions for 
the thermodynamics functions of the same kind found in [37]. On the other hand, in |41j . the 
particle mass is taken as an independent variable. In this case, it is also found a set of consistent 
thermodynamics relations. Other possibility is presented in [28J , where instead of keeping the 
entropy unchanged, as Gorenstein and Yang, the internal energy preserves its SSM definition, 
in a way that the consistency of the thermodynamics relations is also achieved. Finally, a set of 
solutions for the thermodynamics of quasi-particles at \i — was presented in [30J. There it is 
shown that the solutions of references [2E1 E7J are particular solutions of a general formulation. 
The next natural step is to extend this formulation to the case of finite 11. 

In this work, a framework for the thermodynamics of a system composed of quasi-particles 
whose masses are T and /i dependent is developed, generalizing [3D]. At first sight, a system 
composed of a mixture of fermions and bosons should follow the law of Dalton for gas mixtures. 
This law states that the pressure exerted by a mixture of gases is equal to the sum of the partial 
pressures of all the components present in the mixturq^. This statement holds exactly for ideal 
gases. For real gases, however, the equality is approximated. As the quasi-particle approach 
is built on the ideal gas model, one would expect that the pressure for the whole system, i.e. 
for quarks and gluons, were the sum of the fermionic and bosonic parts, separately. In reality, 
this should not be true in the qQGP because the underlining theory is that of an interacting 
system. It makes no sense to fix the number of gluons n g = dP g jdii = 0, since no energy 
is expended for their creation (n g = 0). On the other hand if the pressure is computed from 
Ref. [30], with a suitable Debye mass, mo = m D (T,fi), and including fermions, it can be 
seen from the thermodynamics relations that the number of gluons will be different from zero 
n g = dP g (mr>)/dn ^ 0. Therefore, the law of Dalton fails here, implying that it is necessary to 
develop the theory for the whole system, and not for each part independently. 

This article is organized as follows. In Sec. II, the general formalism for a consistent ther- 
modynamics of quasi-particles in the grand canonical ensemble is presented. In sections III and 



1 This stretch is a quotation from E. Fermi in [42]. 
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IV, particular solutions for this thermodynamics are computed, and in section V a comparison 
with perturbative QCD is made. In section VI, the comparison is with lattice QCD and, finally, 
the conclusions are presented in section VII. 

2 The Formalism 

In the quasi-particle approach for the quark-gluon plasma, the system is described by an effective 
Hamiltonian with mass terms that are dependent on T and/or /i [231 I2H EDI 33j . As widely 
known, early attempts to implement this idea run in a serious problem: an inconsistency in the 
thermodynamics relations. As shown by Gorenstein and Yang [37], it is not enough to use the 
Hamiltonian of an ideal gas with the mass m replaced by the effective mass m(T, /i). An extra 
term, B, has to be added to the original Hamiltonian: H — > H+rjB, where r/B is regarded as the 
zero point energy of the system. The extra term modifies the usual thermodynamics functions, 
healing the thermodynamics inconsistency of the system. Thus, the most general form for the 
thermodynamics functions can be written as: 



$(V,T,fi,m 2 f ,7r%) = $ f (V,T,ii,m}) + MV,T,v,ml)+aB(V,m 2 f ,m 2 b ), (1) 
U(V,T,fi,m 2 f ,m 2 b ) = U f (V, T, /i, m 2 ) + U b {V, T, //, m 2 ) + V B(V, m 2 , m 2 ), (2) 
S(V,T,fi,m 2 f ,m 2 b ) = S f {V, T, fi, m 2 ) + S b (V, T, m 2 ) - ^B{V, m 2 , m 2 ), (3) 

N(V,T,fi,m 2 f ,m 2 b ) = N {f} (V, T, //, m 2 ) - -B{V, m% m 2 b ), (4) 

where $ = —PV is the grand potential, U is the internal energy, S is the entropy, and N is the 
average net quark number. The subscripts / and b denote fermions and bosons, respectively. The 
functions m 2 = mj(T,fi) and ml = m 2 (T,p,) are the dynamical masses acquired by fermions 
and bosons. The constants a, rj, 7 and A are arbitraries, but constrained by a = rj + 7 + A, 
to guarantee the consistency of the fourth thermodynamic relation in Eqs. (jSJ) below. This 
condition is just the extension of the constraint obtained in Ref. |30j. The T and /1 dependency 
of B(V,m 2 f,ml) is assumed to be through the masses only. The extra term B(V,m 2 ,m 2 ) must 
be such that the thermodynamics relations 



'--(£)„• *--(£)„• (5 > 

are valid. Substituting the Eqs. fll}, fl2J), (j3J) and (j4]) in Eqs. (jSJ), one finds the following pair of 
differential equations: 



Using the connection between SSM and thermodynamics, $j = — TlnZj, where Zj is the 
grand partition function with Hamiltonian ifj, one can deduce the following equalities: {^jf) = 
^-(§^)t,v,^ and (^) = ^-{^)t,v^- Using these expressions, and assuming a/0, one 
concludes that the two equations in ([2D are the same. They are rewritten as: 



For a = 0, the limit a — > in Eqs. ([6]) has to be taken with care. One obvious point is that 
the two equations for B in (J6l) give different results if one simply sets a = Co When the grand 
potential has no extra term, the average number of particles and the entropy must have extra 
terms, i.e. if a = then 7 and A must be different from zero. If such condition is not fulfilled 
than the thermodynamics functions will not be consistent. This conclusion is based on the fact 
that the entropy and the average number of particles are obtained from the grand potential by 
Eqs. ([5]). Therefore the previous redefined thermodynamics functions are not valid when a = 0, 
and the modified thermodynamics functions must have the following form: 



$ = $ b + $ f ,S = S b + S f ->y-£,N = N {f} -\-?!-,U = U b + Uf-\B N -> r B s , (8) 
where B$ and B^ are 

2 There is one trivial possibility where B is the same, when 7 = A and = C{F [i. But this is not the kind 
of solution that is expected for the quasi-particle mass. 
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The solutions for B given by Eq. (Q are inversely proportional to A and 7, respectively. Thus, 
Eqs. (IE]) are independent of 7 and A, differently from the a^O case. 

As in the canonical case [30J, the physical meaning of B(V,T,p) can be extracted from 
quantum statistical mechanics. Using the quantum form for the internal energy, one has 

U=(£H)+r,B = Tr[p(J2 H + 17% (10) 

bj b,f 

where p = e~ f3H+f3flN /Z is the ensemble density operator, H is the Hamiltonian operator, N is 
the number operator, and B = Bl is the extra term with the unity operator. Note that when a 
term independent of q and p, but dependent on T/p, is added to the Hamiltonian, the density 
operator for this modified Hamiltonian is exactly the same as that obtained using the original 
Hamiltonian, i.e. the density operator is defined up to a T/p function in the Hamiltonian. Thus, 
one can add B to H in p without inducing any change in Eq. (IT01) . Using this property, it is 
seen that Eq. ffTUj) has the same definition as in the SSM case: U is given by the average value 
of the system Hamiltonian. Therefore, the complete Hamiltonian, particles plus vacuum, which 
must be used in the quasi-particle picture, is Ht = J2H + E , with E = r/B the zero point 
energy. The standard procedure of discarding the zero point energy can not be used here, since 
it results in uncompleted thermodynamics functions, causing an inconsistency. 

For the average number of particles, the same manipulation (as for the internal energy) 
can be applied to the iV quantum expression, where instead of adding a T/p function to the 
Hamiltonian, one adds it to the N operator. Thus, N follows the same definition as in the 
SSM: N is given by the mean value of the total number operator, where the total number 
operator is Nt = N + No, with iVo = —j^B. The extra term —XB/p represents the number of 
particles associated with the zero point energy, which means that one can regard the zero point 
energy of the quasi-particle scenario, multiplied by suitable constants, as a mechanism to add 
or subtract particles from the system. 



These redefinitions imply the following thermodynamics functions: 



S=-(lnp>- 7 _ <$> = -T\nZ T + 1 B, N = (N T ), U = (H T ). (11) 

where Zt is the partition function with Ht and Nt- Note that for 7 = 0, the statistical 
mechanics definitions for the entropy, the grand potential, the number of particles and the 
internal energy are the usual ones. 

The interpretation for B in the a = case follows closely that of the a ^ case. Using 
again the quantum version of Eq. (jSJ), the internal energy can be written as 

U=U b + U f -\B N - 1 B S = Tr[p( £ R { - XB N - 7 B S )], (12) 

i=b,f 

where Bj = Bjl is the operator for j = N,S. As before, Ht = Y^H + Eq, where Eq = 
—XB N — ^Bs is the zero point energy, dependent on T and \i dependent. In the quantum 
statistical mechanics picture, the self-consistent definition for the number of particles is 



iV = (N) — -B N = Tr 
f 1 



p[N- -B N 



(13) 



Applying the same arguments as before, the total number of particles operator is Nt — N + No, 
with A^o = —j^Bn the number of particles associated with the zero point energy. Finally, the 
quantum statistical mechanics relations are written as a function of the total Hamiltonian 

$ = -T\nZ T + jB s , S = -(lnp>-7^, N — (Nt), U = (H T ). (14) 

To conclude this section, the solutions for both a^O and a = cases are summarized in the 
following table: 

3 Thermodynamics Functions for a^O 

To investigate the possible solutions for the thermodynamics of the system, one needs an ex- 
pression for the extra term B. From Eq. (j7j), it is written as: 
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a.B = B jj,^T^ — /j<*Ta ^2 [dm?/ — — «T «, 

i=6,/ ^ \ / 

or simply 

a5 = Bofi^T^ — /i"T« E 

»=&,/ 

As the approach developed here is completely general, particular choices for the constants a, 
7, A, and r], should reproduce known results in the literature, as for instance the equations 
obtained by Gorenstein-Yang [37], and the ones obtained by Bannur [33] . 

3.1 The Solution of Gorenstein-Yang 

In Ref. [37J, the thermodynamics of a system whose masses depend on T j \x was studied through 
two conditions: {%) the thermodynamics functions are defined as in the standard Statistical Me- 
chanics, and {%%) the pressure is minimized with respect to the phenomenological parameters. In 
the general solution presented here, the condition 7 = guarantees (i), implying that the vac- 
uum does not contribute to the total entropy of the system. The corresponding thermodynamics 
functions are then given by: 



, dHA a_ 7 r^/dHA a_ 7 
dfil -t—^ )fT«T~ + dT{ )n"T~* 



dfx 



dT 



(15) 



S 
$ 
N 
U 



i=b,f 

- \ T\nZi + fi«J d/ji 

8=6,/ \ 



a 



E 

i=b,f L 



i=b,f 
V ( A 



\ di* 

IdHi 



)H 



[H^ Lu« / dfj, 



a 



dT I 
IdHj 



+ fj," / djjt, 



)H~ + / dT 



/dHA x 
dH\\ 



aB ^ 



dT 



B \fi<* 



+ rjB fj,c. 



^-1 



(16) 



The solutions with 7 = will be classified as of Gorenstein-Yang type of solution. They can 
be divided into two different branches. 
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3.1.1 The Original Solution 

This is the most used solution, which satisfies both conditions (i) and (ii), and was used in 
several studies of QGP by quasi-particle model, as for instance in Refs. [23J, [2U [251 12SI 127] . 
In this approach, neither the entropy nor the average number of particles are modified by any 
additional term coming from the vacuum. This solution is obtained by setting 7 = A = 0, and 
a — r\. The thermodynamics functions are then reduced to: 



s = - E < ln A> 

i=b,f 



N = (N) 




As this solution is well know and has been largely studied, we will not provide here further 
details of its features, which can be found in Ref. [37]. It will be referred from now on as GY1. 

3.1.2 The Modified Solution 

A second solution of the type Gorenstein-Yang, is the one where a = A and 7 = 77 = 0, meaning 
that both the entropy and the internal energy retain their original forms, while the other 2 
thermodynamics functions are modified. Including bosons and fermions, the thermodynamics 
functions are now written as: 



S = - £ (In ^) 

i=b,f 

u = E*> 

i=b,f 

10 



This solution will be referred as GY2. 
3.2 The Solution for 77 = A = 

It is based on the principle that thermodynamics quantities that have microscopic analogue 
must be computed by the average value of its c/ass«c[] counterpart. For example, the energy 
and the average number of particles. Bannur studied this case for \x = in Ref. [28], and recently 
extended it to /i ^ [33] • In both cases, it was found good agreement with lattice data, except 
in the region T c < T < 1.2T C . Setting 77 = A = in Eqs. (fH|) and (|T5l) . one obtains: 



N = (N) 

u = E(^) 

i=b,f 




3.3 Physical Constraint on B 

The extra term B, sometimes called bag energy or bag pressure, plays the main role in the 
consistency of the thermodynamics relations. Besides this, B must also be physically acceptable. 
The theory must guarantee that any physical solution containing B must recover the physics 
already known, as for instance reach the Stefan-Boltzmann limit at high temperatures, or the 
perturbative results of QCD at finite temperature and/or chemical potential. In particular, in 
the limits where the masses are independent of T and /x, B must be, at most, a constant: 

lim B{m g ,m f ,V) = B . (20) 
m g — ► 
m q -> to 

Using this condition in expression (|15|) . one obtains: 
3 The word classic here means without vacuum contributions. 
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lim B(m g ,m f ,V) = a~ l B T«^. (21) 
m g — ► 
m g -> m 

From this result it follows that Bo must be zero if 7 7^ and A 7^ 0. 

Usually, B is calculated from a differential equation deduced from the Maxwell relation, 
dS/dfi = ON /dT. This procedure leads to an equation in the running coupling and in the extra 
term 5 (23 E21 ES] : 



«t(T, n,g)— + ctfj,{T, (jl, 9 2 )-q^ = b(T, //, g). 

Thus, the knowledge of the coupling and B at \l = 0, and at an arbitrary T, is enough to 
determine these functions at any T and fi. However, as done in this work, B is given by the 
solution of Eqs. (jSJ). Naturally, one should expect both results to give the same answer. In 
fact, they give. The general solution presented here was built to be consistent. If one uses the 
full expressions for the thermodynamics functions to verify whether the Maxwell relation is 
respected or not, one finds that they indeed respect it, without any additional condition on the 
masses or on the coupling. Hence, one can do as in [33] , and borrow the masses of the particles 
(and the coupling) from any given framework, without spoiling the consistency of the theory. 
The choice of a given mass for the quasi-particles will induce a particular form for B, leaving 
the whole theory self-consistent. However, as will be shown here, this consistency holds only for 
the full functions, and not in a perturbative expansion, with one exception only. 

3.4 The Asymptotic Behavior of the Solutions 

The thermodynamics functions will be computed for a system composed by quarks, anti-quarks 
and gluons, with a Hamiltonian H = tUi = \Jk 2 + m 2 (T, fx), where \k\ is the magnitude of 
the quasi-particle momentum vector, and mj(T, //) is the quasi-particle mass. To obtain the 
asymptotic behavior of the solutions, it is necessary to know the quasi-particle mass as a 
function of T and //. As statistical mechanics does not provided this dependency, the standard 
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procedure using the HTL approach for the computation of quasi-particles masses will be 
adopted here. In next-to-leading order (NLO), the mass for gluons and quarks are [2H 1221 [21] : 



m 



m„ 



D 



71 



V2 



g 2 Tm g , 



8N C 



2 / ,,2~ 

T 2 + 



N g9" lrr.2 , 



7T^ 



(22) 



where N g = (N 2 — 1), N c is the number of colors, iV/ is the number of flavors, and m| 



D 

^ 2Nc+N f T 2 + ^^2) is the Debye mass [2H [221 EH]. Including the bare quark mass, m , one 
writes the quark mass as m q 



mg + fn^. For the numerical calculations presented in this 
work, the value m /T = 0.4 will be used, according to Ref. [5]. Using the approximations g<Cl 
(and n/T <^ 1), the quark mass is then written as 



8N C 




3N f 



N c g 2 



7T S 



(23) 



6 y/l2N c + 6N f * 2 T 2 t 

The coupling constant g will be replaced by an effective coupling constant g s [32], inspired by 
perturbative QCD, in order to accommodate non-perturbative effects at the vicinity of T c when 
fitting the lattice data. In that case, 



16tt 2 



20! ln(ln£ 2 ) 



(24) 



2 m z 

T3 



\ Is. 



with /3 = (UN c -2N f )/3, fa = (3AN 2 - 13N f N c + 3N f /N c )/6, and f = A s £ ^1 + 
where A s is a scale parameter and T s is a temperature shift that regulates the infrared divergence 



below the critical temperature T c . The factor y 1 + u 2 ^ i n £ comes from the study of QCD 
running coupling at finite temperature and quark chemical potential based on semiclassical 
background field method [43J. The parameter u will be chosen to provide a best fit for the 
lattice data. For /i/T — > 0, Eq. fl24"|) is reduced to the usual result for the coupling constant. 

The solutions for a ^ are not easy to calculate explicitly since it is necessary to perform 
integrals on both T and \i to find B. Nevertheless, in the limits T — > 00 and /i — > 00 calculations 
can be done. Eq. (|T5j) is rewritten as: 
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= B b + B f , (25) 

where B b is the piece correspondent to gluons, and Bf contains B and the quark part. The 
average value of the derivatives of the Hamiltonian (with respect to the masses) are computed, 
and the results are: 



D Bq X ; 

a 



a 



E 

b.f 



- T - 




~'T~ 



dH b \ Vd b f k 2 + d b VT 2 (m b \ d b VT 2 

dk , nt(k) = -^-Js -=r ~ — ^j— ! ( 26 ) 



dm 2 / An* J .K^tf h b "' 2tt 2 *\Tj 24 



/dHA _ Vd f r k 2 ^ _ d f VT 2 jT (mf ^ ^ d f VT 2 (n 2 £ \ 



where d b and df are the degeneracy factors of bosons and fermions, respectively. The functions 
If and Jj are given in Appendix A. The derivatives of the mass with respect to T and [i are 



dm 2 2 N f p dm 2 2 2N C + N f 

~W = 9 T^' ~dT = 9 6 T ' (28) 

9// ~ 5 4A c vr 2 ' <9T ~ P 4A C 1 yj 

Note that the asymptotic approximation implies that the masses are calculated in leading order. 
With these four equations it is easy to compute B. For the gluon contribution, using Eqs. (125]) . 
and (EHD, one has 



Btf* „ rn 2 T 2 2 (2N C + N f )T* 2 + A - V 

Vd h ~ 12(2 — A) 12 2 (2-A)(4-T)' K ' 

where A = ^ and r = ^. From Eqs. ( 1251) . ( 1271) and ( |29i) . the quark and anti-quark contribution 
is written as 



Bja r 
Vd f ' 



'Bqh T 



m. 



12(4 -r) 



T 



71- 



g 2 Ny 



32tt 2 (4 - T)A C 



A-T ft 2 
4- Att 2 
2 - 



+ 



2-r 



+ 



r + A V 



3(2- A) 



(31) 
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In the calculation of B^ and Bf at high T and /i, the following asymptotic conditions were 
used: ^ < 1, and — C 1, where To and llq are the lower limits of the integrals. Also, in the 
region g — ► 0, g is a very slowly function of T and fi, what justifies somewhat to regard g in the 
mass derivatives as constants, Eqs. ([28]) and (1291) . The integrals converge in this region only if 
conditions A < 2 and T < 2 are satisfied. Inspired by the case of a system made of bosons only 
at \i = [5U] , this condition will be called of weak physical condition. These results provide the 
necessary quantities to describe the thermodynamics of the system. To write the asymptotic 
pressure one needs Eqs. ([T]), (I3"U]) and (13T|) . Then, 



90 1 2 c / vr 2 2 tH 



^ + av, 2 - r y4 + rA-r-A ^ + Nf lzAt\ . (32) 



32 V 9 4-T J (2-r)(2-A) tt 2 M-Att 4 , 

The term independent of g is the Stefan-Boltzmann pressure, while the rest of the expression is 
the correction to order g 2 . Using Eq. (132]) in the thermodynamics relations Eq. (jSJ) is possible 
to obtain the other expressions. On the other hand, following the same steps used to calculate 
the pressure, one can do the calculation of the other quantities directly. As the theory has 
thermodynamic consistency, both ways should give the same result. However, consistency has 
been proved only for the full thermodynamics functions. One question then arises: Are the 
perturbative thermodynamics functions also consistent order by order? To answer this question 
is necessary to compute S, N and U in both manners and compare the results. 

The asymptotic entropy density, s = S/V, is calculated from Eqs. (jSJ), (1311 and (|3~T]) : 



s{T, fj) 



tt 2 T (7N c N f + AN g 



15 



T 2 + N c N f ^\ + 



24 



4/V c + 5iV)2-r 2 , 3N/ 4 - 6r - 2A + 2IA + T 2 /i 2 3A^T fi 4 



4-r 



(33) 



2 (2-T)(2-A) vr 2 2(4-A)tt 4 T 2 _ 

For the asymptotic particle number density, n = N/V, one writes it in terms of the integral 
involving fermions, which can be found in Appendix|A] (N) = f) . The net number 

density is then written as: 
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n(T, fi) 

Me ( n ,e 

16tt 2 \ t 



7T 



2T - 6 A + 2T A + A 2 2 2A^(2 - A) /z 2 4A^ C + 5 Avr 2 T 4 N 



(2-T)(2-A) 



(4 -A) TT 2 9(4 -T) 



(34) 



On the other hand, computing the entropy and the average number of particles from the result 
for the pressure, Eq. (I32l) . using the thermodynamics relations Eq. (J5j), one obtains: 



Sthe 



dP 
df 



n 2 T f7N c N f + AN, 
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'"-T 3 + N c N f ^) + 



n t he 



dP_ _ N c N f fi 



T2 + ^. 



4A^ C + 5Nf 2 — T 



9 



4- r 



T 2 + Nf 



N 9 N f g 2 ^ 



rA-r 



(2-T)(2-A) 



TA-T-A ii 2 ' 

(2 - r)(2 - A) ^_ 

A - 2-Au 2 
T H — 

4 — A 7T 2 



(35) 



3 \~ ' 7T 2 J 87T 2 

Trivially, the resulting calculation for the independent quasi-particle terms, i.e. the terms 
independent of g, are consistent: They result the same expression for both procedures. However, 
the terms in g 2 fails the consistency check. Take, for instance, the expressions for the entropy. 
Only the terms that depend solely on the temperature match when using both procedures. 
The same happens for the particle number density, only now instead of the terms involving the 
temperature, are the terms involving the chemical potential the ones that matches. The terms 
involving both the temperature and the chemical potential do not match in any situation. 
Comparing both procedures for the terms in \xT 2 of n and in T/i 2 of s, the unique solution 
which produces a possible agreement between the two procedures is the one with T = A = 2. 
But this solution was excluded in the B calculation, the weak physical condition, because the 
integrals do not converge. Concluding, the asymptotic limit of thermodynamics functions breaks 
thermodynamics consistency order by order in a perturbative expansion. Does the same happen 
with the a = solutions? 



4 The a = solution 



The main feature of the a = solution is the fact that the extra terms are easier to compute, 
as they do not involve integrals on T and/or \i. In such case, it is possible to have an analyt- 
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ical result for almost all range of temperature and chemical potential. In order to obtain the 
thermodynamics of the system, it is convenient to start with the pressure, as there are no extra 
terms contributing to it. The pressure of an ideal gas with modified masses is 



P(T, n) = — / dkk A ]T d 

D7T JO ~T, 



nt(k) 



tf \/k 2 + m 2 (T, fi)' 



(36) 



where the distribution functions are n^{k) = [e^^ k +rn 'f ^ + 1] 1 ± [e^v^ +m } + ^ + 1] 1 and 
nf(k) = \{[e p W k2+m D - l]- 1 ± [ e PW k2+m D - l]- 1 ). Expression flU can be computed explicitly, 



( l) ^, -^r < 7T Or (11 J 



a z i 2ai 



< 1, 



see Appendix |Aj using convergent boundary for fermions:(i) — — 
where a = £ and 77 = For bosons, the boundary condition is r b = ^ < 2n. The pressure 
can then be written as P = ^T 4 [dfl^(rf, a) + d b J 5 (r b )], or 



vr 2 T 4 



90 



4 + 



7d f 15a 2 d f 15a*d f 15(d f r 2 f + d b r 2 b ) A5a 2 r 2 f d f 15r fc 3 d, 
~4~ + 2vr 2 + 4vr 4 4^ 4^ + 2tt 2 



45 



■in" 



7T 



4tt 4 
r td b /3 



+ 



105C(3)d/ r , , 



45(r 2 - a 2 ) 2 ^ 



64tt 2 



21C(3)a 2 - 



31C(5), 



_2\2 



327T 2 

3 00 

I 

k= 



+ 



A5d 00 



E 



64tt 6 

-i) n r(n 



/ 

2 " b 



+ 



a 2 — r? 



(27r) 2 ™r Q)r(n + 3) 
2ai 



45(2n + l) 3 - rtj+fe) 

J n=o 2 jfer(±)r(* + 3) lj(2n + 1)tt] 2 + (2n + 1)tt, 

For the computation of the entropy, Eqs. (jSJ) and ([9]) are necessary. In this case, 



k+2' 



(37) 



s(T,/i) 



1 



2vr 2 T jo 



dkk 2 ^2 di 

i=b,f 



U 2 +(m 2 -f^ 



k 2 + m 2 



7P 



/ 2 



16^(r,,a) + 



1 dm 2 . 



nf(k) -firii (k) 



lUr h a)-^Gl{r s 



\ T 2 IT dT 



T 



d h T* 



7T- 1 



mi 



1 J,h) 



+ 



(3f 



An explicit integration can be done in this expression for s using the integrals from Appendix 
[A1 The same sort of manipulation can be made for the case of the internal energy density, 
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e = U/V. The result is: 



e(T, fj) 



2tt 2 Jo 
if 



A- 2 



6,/ J /c 2 + m 2 



,2 . 2 1 / ^m 2 . rpdmV 
k + mj a n + 1 — — 

1 2 V dfji dT 



dinf(k) 



d f T 4 



12/ 5 T (r /; a) + 



mj_ 

J^2 



dm 2 . 



2T 2 



+ T- 



_dm 2 
df 



d b T 4 



71^ 



l2J 5 (r b ) + 



mi 



T 2 2T 2 



I'- 



dm\ 



,dm\ 
~dT 



Similarly, the result for the average number of particles is 



(39) 



n(T, n) 



2vr 

IT 2 



I poo 

— / dkk 2 V d 



n i 0) - 



dm. 



2Jk 2 + m 2 9fi 



in+(k) 



d f G%(r f , 



1 

2T 



df ~di I ^ Tf ' a) + db ~diT Mrb) 



(40) 



As done in the a^O case, it is necessary investigate the consistency of the expanded thermo- 
dynamics functions. The pressure at leading order can be obtained from Eq. (|37l) . 



7T 2 7N C N P + AN Q . , r u 2 T 2 N f N c a 4 

6 30 } TT 2 2 7T 4 



N q g 2 AN c + hN f . Ar fi 2 T 2 N f // 4 \ 

18 / 7T 2 2 7T 4 / V ^ 



16 V 

The leading order results for s and n are calculated from Eqs. (158]) and (H0]1 . The results are: 



= w r ^_^^ + g). (42) 



7T 2 / »7T^ \ 7T 6 , 

The entropy density and the particle number density can also be calculated directly from 
the thermodynamics relations. As before, we denote them by s t herm and n t herm, respectively. 
Opposite to the a ^ case, the resulting expressions now agree: For a = 0, s = s t herm and 
n = n therm at order g 2 . 
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5 Comparison with Perturbative QCD 



The quasi-particle picture has been introduced to describe the QGP mainly because perturba- 
tive QCD is unable to study the QGP close to the deconfmement region. Perturbative calcula- 
tions, indeed, exhibit a poor convergence except for temperatures as high as > 10 5 T C . In any 
case, perturbative QCD calculations are what the model should reproduce at extremely high 
temperatures and chemical potential. 

The thermodynamics functions, to order g 2 , in perturbative QCD are [21~1 [22]: 



tt 2 ( 7N c N g + 4N v . fj?T 2 N f N e yT 

6 30 1 c 7T 2 2 7T 4 



N g g 2 ( W c + 5N f /i 2 T 2 N f ^ 

S QCD = ^ (!^±i^ T2 + NcN £\ _ ^ (1M T2 + ^ 



nQCD ^(T 2 + ^)-q^(T 2 + ^. (43) 



It is fundamental that the thermodynamics functions calculated here be compared with 
the results from perturbative QCD. First, the solutions for a ^ are compared. As discussed 
before, these solutions have the serious problem of not being consistent order by order. Never- 
theless, the comparison with perturbative QCD is made using expressions (I3"2j) . (13"3l) . and (1541) 
computed directly from their definitions. If one uses, instead, Eqs. (135|) . the expressions for the 
thermodynamics functions would be different, but the final conclusion would be the same. The 
cases discussed are: 



1. 7 = A = 0, a ^ 0: Solution GY1. 



g 2 N g [4N c +5N frr4 , AT ft 



G 64 V 9 7T y 

g*N g T ( 4N c + 5N f rri2 ^\ g 2 N g N f y ( g 



7T 2 



Although this is the solution widely used in the literature [231 1211 [25], [261 EZ] , the asymp- 
totic expression for the pressure does not agree with perturbative QCD, although s and 
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n agree. However, this agreement is purely accidental given the lack of consistency of 
the thermodynamics functions at this order, as discussed in section 3.4. Note that the 
whole problem arises in the mixed term, fi 2 T 2 , in the pressure. This is the first explicit 
calculation of these 3 functions. 

2. 7 = r] = 0, a ^ 0: Solution GY2. 

g 2 N g (4N c + 5N frr4 N f /i 2 T 2 N f ft 



Gy2 32 I 18 2 vr 2 3 vr 4 



2 2iV 9 T (4N c + 5N frr2 AT f 



sgy2 = -Z-^-l—^; '-T* + N 



16 V 9 ' ' "V J ' 
ff 2 iy ( AN C + SiV/ 7r 2 r 4 N f 2 2N f » 2 \ 
Gy2 16tt 2 ^36 /x 2 2 3 vr 2 / 

Here only s match with perturbative QCD. 

3. A = rj = 0, a ^ 0: Solution proposed by Bannur [281 [33] . 

= ^ + 5^ _ ^ JV^ 

32 V 27 ' vr 2 2 vr 4 / ' 

g*N g T ( 4N e + W f ^ 3N f 3N f // \ g 2 N g N f » ( 2 /i 2 \ 

B " 24 ^ 9 4 vr 2 8 vr 4 T 2 y ' S 16vr 2 ^ vr 2 ,/ " 

As expected, the only function whose asymptotic form agrees with QCD is n. 

4. a = 0: The consistent solution. 

° 16 \ 18 ' 7T 2 2 Vr 4 / ' 

g*TN g ( 4N c + 5N f ^ a 2 \ g 2 N 9 N f » ( 2 ^\ 
Sq=0 = 8~ i 9 T+iV/ ^J' ^= = 8T 2- l v T+ ^J- 

Unlike the 3 precedent cases, the solutions with a = are the on/?/ solutions which are 
consistent order by order. By extension, they are the only solutions that are meaningful 
when expanding in g 2 . However, when using the HTL/HDL masses for the explicit cal- 
culation of the thermodynamics functions, Eqs. (j4ip and (|4"2]) . one gets half the result of 
perturbative QCD for the 3 functions. 
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The relevant question now is: Is there any scheme that is consistent order by order in the 
asymptotic limit and that also reproduces perturbative QCD? As seen, the only possibility of 
consistent schemes at finite T and \x are the ones with a = 0. However, to match perturbative 
QCD, a quasi-mass for quarks and gluons that is not the HTL/HDL mass has to be used. 
Inspired by the calculation at finite T and fi = of Ref . [30] , where the HTL mass was replaced 
by ttihtl/Vz in order to match perturbative QCD in the a = 0, the same can be done here. 
In fact, doing this replacement for both quark and gluon quasi-masses, the agreement with 
perturbative QCD is achieved. 

The next natural question is: Does the scheme with a = also reproduces lattice QCD data 
in the region near T c ? This question is answered in the next section. 



6 Lattice data and qQGP 

In this section is discussed one of the most important features of the quasi-particle approach: 
The fact that it is a useful tool to describe lattice QCD data. Lattice QCD at non-vanishing 
chemical potential uses a Taylor expansion in (/x/T) around // = to compute observable 
quantities up to 0(fi 6 ) [5\. For the case of the pressure, the expression is: 



P 



co(T) + c 2 (T)^ + c 4 (T)^ + c 6 (T)^ + ■ 



(44) 



The coefficients c 2 ,C4,C6 have been computed in Ref. [5] in the range [0.76,1.98] for T/T c . In 
the present approach, the parameter c 2 will be used to constrain the parameters u, X s , and 
T s appearing the running coupling at finite T and /i, Eq. (|2"4"I) . As seen in the last section, 



the only case that reproduces perturbative QCD at is for a = 0. In that case, p 
-\[dflj(rf, a) + o?bJ 5 (rfe)]. On the other hand, the definition of c 2 is c 2 (T) 
the derivatives, one gets: 



J^4 



d 2 (P/T 4 ) 
da 2 



a=0 



Doing 



o 2 p 



da? 



a=0 



d \w a 



d 2 p 



da? 



1 




a=0 Ti-2 





StT 



3 da? 



d b r b Ji 



d 2 r b 
da 2 



a=0 



a=0 



(45) 
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The expressions for If, /J\ and J3 have sums that have to be stopped somewhere when doing 
the numerical calculation. The criterion used was the radius of convergence of these sums, which 
is r b j < n. With this in mind, one sees that for r b j close to ir, orders bigger that k = 5 are 
negligible. 

The data points from Ref. [5] where corrected by a factor of 1.12, because as stated in Ref. 
[46J, there should be a correction of 10-20% when extrapolating the data to the continuum. 
Following Letessier and Rafelski [47J, the first choice for the parameter u appearing in the 
coupling was u = 1/tt. However, with this value, the fit for the data below T/T c = 1.2 is 
quite poor. If one chooses, on the other hand, u = l/2ir, agreement between lattice and the 
calculation presented here is satisfactory, as long as 2.7 < A s < 3 and 1.6 < T s < 2. With these 
results the reduced pressure = P( - T ' / ^ 4 P( - T ' - > can be calculated and compared to the data. 
This is done in Fig. 1, with X s = 3 and T s = 1.98. As seen, the agreement is quite good. 

7 Conclusion 

The number of different approaches that results in a consistent thermodynamics of quasi- 
particles, whose masses are T and fi dependent, is large, as the general formalism presented 
here shows. The question is which one of the possible solutions is the most adequate. Or, in 
other words, which one allows one to get the most from it with some minimum input. Clearly, 
one can say that it is an advantage to work with simpler expressions, preferably with equations 
that allows one to make algebraic manipulations all over the calculation. And, of course, the 
adequate solution should be the one that reproduce both perturbative QCD and lattice QCD. 
Besides, it is desirable to have a solution with a field theory analogue, as for instance the $ 
derivable approximation, in order to frame it in a stronger theoretical level. 

From the solutions studied here, a whole class of them (a ^ 0), which includes the Gorenstein- 
Yang type of solutions, have problems when a perturbative expansion is made. Specifically, al- 
though the thermodynamics consistency is achieved for the full solution, it is lost order-by-order 
in the perturbative expansion. This implies that this class of solutions are not able to reproduce 
perturbative QCD at finite T and \i in any scenario. Of course, when trying to reproduce lattice 
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data near T c , such problem does not appear because, in that case, the full solution from the flow 
equation for B is used, and no order-by-order approximation is made. Despite these problems, 
the solution GY1 has on its side the fact that the entropy density and the particle number den- 
sity preserve the same form of the standard statistical mechanics, exactly as in the $ derivable 
approximation. On the other hand, for the class of solutions with a = (pressure is unchanged 
by the extra term B while the other functions are modified), order-by-order thermodynamics 
consistency in the asymptotic limit is maintained and agreement with perturbative QCD is 
achieved. Moreover, the lattice QCD data is reproduced down to T/T c ~ 1.02. 

In summary, a general consistent approach for the thermodynamics of quasi-particles at 
finite T and \x was presented for the first time. It was show that the only class of solutions for 
the thermodynamics functions capable to reproduce both perturbative and lattice QCD is the 
one where the pressure is not modified by the introduction of the extra term B. 

This work was supported by FAPESP (04/15276-2) and CNPq (307284/2006-9). 

A Appendix 

The necessary integrals that appear throughout this work are: 

A.l Bose-Einstein Integrals 

The Bose-Einstein integrals are [3U1 HE] 



X 

J n (r) = ~ , / dx 5 , n > 0. 

T(n)J (a? + r 2 )3 e (*»+r»)* _i 

The relevant integrals for our calculation, following [HI H5j [30] are: 



_ , tt 7r r / r 1\ 1 ~ . 

Jstr) = r In h 7r > (m + 2)r 

dV ' 12 4 8 V 4vr ' 2 J 4^ V u 



r 2(m+l) 
m' 



^ 2* 2 * 3 ^ A r 3\ 1 ~ . 2(m+2) 



2 3 425 96 48 2 7 V 4tt ' 4/ 32 m , 
where j m = 22m^^2mr(Y/ 2 2 \r(r+2) • Tne series are convergent for r < 2n. 
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A. 2 Fermi-Dirac Integrals 

The Fermi-Dirac integrals are: 



I roc X U ^ 1 

I„(r,±a) = ——- / dx . = — , „ = , n>0 

V ; r(n)Jo Vx 2 + r 2 e ^+^ ±a + l 

I poo X n ^ 

G n (r,±a) = —. . / dx — , = = . n>0 

v ' ; r(n) Jo e V^+^±a + ! 

As the system is composed by quarks and anti-quarks, the thermodynamics functions 
written as a linear combination of these integrals. Thus, the total integrals are: 



In(r,a) = I n (r,a) + I n (r,-a), 
G l( r , a ) = G n (r,-a)-G n (r,a). 

These integrals satisfy the following recursion relations: 



dlj l+l (r,a) 

dr n 
dl% +1 (r,a) _ 1 T 



--/J_i(r,a), 



da ~ G n-^a), 

or n 
<)(■!,. t (r. a) 



nl n+1 (r,a) + —I^faa), 



da 

with the initial conditions 



2(1 — 2 2 ~ n ) 

/J(0) = -i J -((n - 1) n > 2, and G£(0) = n > 



n- 1 

where is 



m=0 ^ 
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the Riemann zeta function. The procedure to compute the integrals is inspired by Ref. 
where bosonic integrals in the high temperature limit were computed. 

If the integrals If and Gj are calculated, the recursion relations can be used to determine 
the remaining integrals. For the integral Ij , the following identity is used [48J : 



tanh(^) = 2 £ 



to z 2 + [(n + \ )7rp 



1 - 



e 22 + 1 



Thus, the Fermi-Dirac distribution is written as: 



1 1 ^ 

2 ^t2j_ IY9« J- 1 W12 ' 



- x 2 + [(2n + 1)tt] ; 
Substituting this expression in J^, with x — > (x 2 + r 2 )^ + a, follows 



J 1 r (r,a)= / 



Vx^T 



-2E 



Va; 2 + r 2 + 



VrT? — a 



V (Vx 2 + r 2 + a) 2 + [(2n + l)vr] 2 (Vx 2 + r 2 - a) 2 + [(2n + 1)tt] 2 , 
Manipulating the integrand to simplify the integral, one gets: 



if(r,o) 



dx 



n=0 



x 2 + r 2 — (a — (2n + l)7ri) 2 



ii-4^i 2 , n . (53) 

n=0 



To solve the sub-integrals i\ and Z2,n, it is necessary to know the solution of two specifics 
integrals: 



dx- 



x 



o ~~ Vx 2 + 6 2 2 V 2 '2/' 



(54) 



w 



J 

Jo 



OO r£ 

dx 



x 2 + 6 2 26 x + 



1 fe + l 1-e 



(55) 
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where B(p,q) is the Beta function. The factor x 6 was introduced as a convergence factor. At 
the end of the calculation, the limit e — > is taken. Hence, the sub-integrals are: 



OO r£ 

dx- 



o Vx r + 



r " /l — e e 
2 2 



(56) 



z 2 ,n = / dxR — — . ,, 2 

jo \x z + r z — (a — {2n + l)irt) z ^ 



Thus, the integral lf(r,a), Eq. (1531) . is given by: 



ft 



° \ 2 1 : 



2[r 2 -(a-(2n + l)7ri) 2 ] 2 



(57) 



" r V ' ' 



If (r, a) 



i 



2 V 2 '27 sin[f(l-e)]^ \[ r * - (a + pn+l)™)*] 1 ? 
For e — > 0, the functions can be expanded in power series: 



(5f 



B 



b~ € = e 
1 - e e 
2 '2 



In 6 



CSC 



2^ ; 



1 -eln6 + C(e 2 ), 

^ + In 4 + f y + In 2 2^J e + £>(e 2 

l + 0(e 2 ). 



(59) 



On the other hand, the denominator in the sum is rewritten as 



1 1 

[r 2 - (a + (2n + l)™) 2 ] 1 ^ ~ [(2ra + l)^ 

1 ~ T(^ + fc) / a 2 2ai r 2 \ fc 

+ [(2n + 1)7^ ^ r(i±£)T(A; + 1) \ [(2n + 1)tt] 2 + (2n + 1)tt [(2n + 1)tt] 2 J ' 

where a Taylor expansion was done. This expansion is valid for two cases: a, r < ir and 

I ( [(2n+l)^ + T>ST> - [(2n+l)H 2 ) ' < L Umit C ^ ' OIle getS 

00 1 2 1+e — 1 1/1 2\ 

£ [(2n + IKF = T^ C(1 + C) = ^ U + 7 " + ^ n) ' 
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where 7^ m 0.5772 is the Euler-Mascheroni constant. The final result for If is then 



11 ; t 7 ^2n + i^r(|)r(* + i) 



a 2 — r 2 



2ai 



[(2n + l)7r] 2 (2n + l)7r 



(60) 



The integral If was also computed in [43J, using the Bose integrals of Ref. [IB]. In [13], however, 
the integral is truncated at some order. In the present calculation, there is no truncation. 
To compute the Jj(r, a) integral, one needs Gf(r,a). Using Eq. (ji7jl 



\Jx 2 + r 2 + 



Crf(r, a) = Gi(r, — a) — Gi(r, a) 

/■oo 00 

= 21 dx V" ^=^_ 

io ^ [(VPT^+a) 2 + [(2n + l)7r] 2 

After some manipulation, 



r" — a 



( Vx 2 + r 2 - a) 2 + [(2n + 1)tt] 2 . 



(r, a) 



n=0 



a — (2n + l)7ri 



x 2 + r 2 — (a — (2n + l)m) 



-4 E 

n=0 



The «3 jn integral is similar to i2, n - 



(61) 



3? 



oT e (a - (2n + 



B l+e 1-e 
1 2 ' 2 



x 2 + r 2 — (a — (2n + l)7u) 2 
As before, taking the limit e — > 0, one gets 



3? 



a — (2n + l)7rz 



[r 2 -(a-(2n + l)7ri) 2 ]^ 



GT(r,a)=a-^^-^C(3) 



2 (27T) 1 



a 2 — r 2 



2ai 



■ a + (2n + 1)ttc " T(| + fc) 
n 4f 1 2n + l fer(i)r(fc + l) l v [(2n + l)7r] 2 + (2n + l) 7 r i 



(62) 



Notice that it was necessary to use both /c = and = 1 terms in the series expansion to get 
a finite integral. Using the recursion relations (|48[) and (|49p . with Eqs. (I6D1) and (|62p . one can 
easily integrate both expressions and use the initial condition Eq. (I52~I) to obtain Ij: 
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rp , . 7i cb r / 1 r\ 7(r 2 -a 2 ) 2 , 

/ T r, a = — + - 7s - In ^- — ^( 3 + 

3V ' ; 12 4 4 V2 ' n J 16 (2vr) 2 sv ; 

_ - 7r 2 (2n + l) « r(| + fc) / a 2 -r 2 2az \ 

n ^ 2 ^r(|)r(fc + 2) U(2" + 1)t] 2 (2n + l)7rJ " 1 1 

To compute Gg one needs Eqs. (150]) and (I5T1) . One starts taking the derivative of the solution 
of Eq. fl50|) with respect to a. Then after some manipulation of this result and comparing it to 
Eq. (lo"Tj) . one gets 



7T 2 a a 3 



6 f 6 



ar 



7C(3)a 
16(2tt) 2 



• 2 2\2 i 

r — a ) + 



21C(3)a 



16(2tt) 



^[(2n+ 1)tt] 2 K 

n=0 



( 4a V - r 4 - 3a 4 



31C(5) 
32(2vr) 4 

r(| + fc) 



2\3 



/ 2 2\ 

a(r — a ) 



+ 



a 2 — r 2 



2ai 



a + (2n + 1)tu ^ 

2(2n + 1) ^ 3 r(i)r(fc + 2) V [{2n + 1)tt] 2 ' (2n + l)vr < 



k+l' 



.(64) 



The second line is the contribution from k = 2. Once Gf 1 is known, can be determined. From 
Eqs. (HHD and ([MD, and from Eqs. (|39|) and (jMjh one has: 



7 5 J (r, a) 



7tt 4 
2 5 45 



7r 2 (r 2 - 2a 2 ) a*_ 
96 + 96 



a 2 r 2 



r 

32 + 43 



3 7" 
7 - 1e - In - 

4 7T. 



7C(3) 
2 7 3(2tt) 



a 2 ) 3 + 



+ 



21C(3) 2 . 2 2 . 2 31C(5) . 2 



2 7 (2tt) 



2/ 2 2\2 

-a (r — a j — 



2 10 (2vr) 4 



r — a 



2\4 



+ 



~ vr 4 (2n + l) 3 » r(| + A;) , |X 

L — is — ^ - ' h 



a — r 



ra=0 



+ 



2 cm 



fc+2 



fc=2 V2 



r(i)r(fc + 3) V[( 2 ^ + l)vr] 2 (2n + l)7r 
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Figure 1: The reduced pressure Ap/T 4 . The points have been obtained from lattice QCD [5]. 
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Table 1: Summary of the two kinds of solutions for thermodynamics relations. 
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